      function chstar(xmu)

      implicit real*8(a-h,o-z)
      chstar=3.2165+4.3431*xmu+.536*xmu*xmu-.0781*xmu*xmu*xmu
      return
      end